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Abstract 

Distribution function and current density in a one-dimensional super- 
lattice with parabolic miniband are calculated. The current dependence 
on the temperature coincides with experimental data. Generalization is 
carried out to quasi-two-dimensional superlattice with paraboloidal mini- 
band. For a sample opened in Y direction with dc current in X di- 
rection, a novel nonequilibrium phase transition is found, namely, ap- 
pearing a spontaneous transverse electric field E y under temperature ris- 
ing. Near the transition temperature Tc, determined by the applied field, 
E y ~ ±VT - T c . 

1 Introduction 

In present work, we study spontaneous appearance of a transverse (with respect 
to the current j x in the sample) electric field E y in a quasi-two-dimensional su- 
perlattice (2SL). In fact, a nonequilibrium second-order phase transition (NPT2) 
is considered, where E y takes a part of an order parameter, while driving field E x 
is a control parameter. Such a phenomenon is similar to the multi- valued Sasaki 
effect in multivalley semiconductors, that was observed in experiments .1 . That 
effect, however, has been not treated as a NPT2, meanwhile such an approach is 
fruitful enough. This approach allows to develop fluctuation theory of NPT and 
reveal noise-induced NPT, as well to investigate stochastic and vibration reso- 
nances in 2SL and other systems Unlike the earlier published works [2]- 
jS], where the conventional cosine- type model was used for the conduction mini- 
band, here a dispersion law is considered in form of a truncated parabola, i. e. 
the dispersion law is assumed to be parabolic up to the Brillouin zone edge. 
Such a problem statement is of interest, because the modern technology allows 
to vary widely the form of the potential relief and the SL energy spectrum. 
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In the low temperature limit (T — > 0) the 1SL conductivity with parabolic 
miniband was studied in ,6j . Here we come out the scope of the approximation 
used in |S] and find the temperature dependence of the current density in the 
model mentioned. It appears, that other NPTs, as compared to |2| _ |3], are 
possible in 2SL with parabolic miniband. In those NPTs, the temperature is 
the control parameter, so that a transverse e.m.f. appears spontaneously under 
temperature raising. 



2 Distribution function and current density in 
1SL 

The electron energy in the 1SL lowest miniband is ^6 

P 2 

e(p) =e(p±) + ^-, -nfr/d < p x < nh/d, (1) 
2m 

where p is quasimomcntum, d is SL period, X axis being directed along the 
SL axis, s(p±) is in-plane electron energy, ir 2 h 2 /md 2 = A is double miniband 
width. 

In quasi-classical situation (A 3> eEd, H/t; t being electron momentum 
relaxation time, e is electron charge), the current density in uniform dc electric 
field E is found by solving Boltzmann equation with collision integral within 
r-approximation: 

sE,^ (p) V f ° (p) ~ f(P) , (2) 



dp J 

where F (p) is equilibrium electron distribution function, Ftp) is unknown dis- 
tribution function perturbed due the electric field. 

We use dimensionless variables below by changing irhp/d — > p, E/Eq — > 
E, T/A — > T (Eq — irh/edr, T is temperature in energy units). 

With the field E is directed along the 1SL axis, we have F(p) = fo(p±)f(p x ), 
-Pb(p) == /o(p±)/o(Px)j where fo(p x ) is equilibrium distribution function, nor- 
malized to the carrier density n (/o(p±) being normalized to unity). Thus, 
fijpx) function satisfies the following equation: 

E x ^- = h-f (-Kp B <l). (3) 
dp x 

We consider nondegenerate electron gas, so that 

f (p x ) = 2n L/WTevi (^=)) cxp f-^j (4) 

(erf(z) is error function). 

In the low temperature limit (T — ► 0), Eq. Q reduces to the function used 
in® 

go (Px )= 2n5(p x ). (5) 
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The exact solution of Eq. takes the form 
1 



/(Px) = 



E 



*erf (l/ 



exp 



E x 2El) 



erf 



T 



+C exp 



Px 
E, r 



(6) 



The second term is general solution of the homogeneous equation J2J). The 
constant C is found from the function f(p x ) periodicity, /(— 1) = /(I)- Then 
we get 



/(Px) 



£ x erf 1/V2T 



exp 



/.Pa 



[exp(2/ J B ;r )-l] 'erf 



£ 5 

T 1 



2£?v' 



erf 



/ Px 



•[l-exp(-2/£^)] 1 erf( 



T 



2T 
1 



V2^ K 



2T VIE, 



-1 <p x < 1. 



(7) 



In limiting case E x — > Eq. J7J) reduces to Eq. J3J|. In another limiting case, 
T — > 0, we get the distribution function found in [H]: 



5(Px) = -=-exp - — 



J[l-exp(-2/^)]- 1 , 0<p :c <l ) 



£j \ [exp(2/.E x )- I]" 1 , -K Px <0 



(8) 



The function J7J satisfies the same normalization condition as the equilib- 
rium function fo(p x ) 



\ J f(Px)dp x = n 



(9) 



and, therefore, it makes the collision integral vanish. Besides, the left-hand side 
of the Boltzmann equation @ vanishes too, because of the periodicity condition 
mentioned. The condition @ can be proved by direct calculation using formulae 
from 0. The distribution function f(p x ) at several values of E x and T is shown 
in Fig. Q] 

The current density (in dimensional units) can be found by a conventional 
way: 

nh/ d 

ed 



2nhm 



P x f(p x )dp x . 



(10) 



By substitution function into Eq. (|10H we get 

j x =E X + [2erf (1/V2T) sinh(l/£ x )] _1 exp [^) 



erf 



T 



\V2E X 



s/ztJ V 



T 



1 



V2E X 



(11) 
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Figure 1: Distribution function f{p x ) &t various values of the driving field and 
temperature, a — E = 0.1, T = 0.1; b — E = 1, T = 0.1; c — £7 = 3, T = 2. 



Here ja, is expressed in units of jo = neAd/nh, while all the quantities are 
written in dimensionless form. Equation (|ll|l determines the current-voltage 
curve (CVC) for the parabolic miniband 1SL with the current density temper- 
ature dependence taking into account. CVC behavior at various temperatures 
is shown in Fig. |2] (CVC at T — > was shown in Fig. 1 of [H]). 



O.OS 
0.06 h 




Figure 2: The current density at various temperature values: a - T — 0, 5; b 
T = 1: c- T = 2. 



In limiting case T — > Eq. I|ll|l reduces to the expression found in j^j: 
j =E 1 



(12) 



svak{l/E x ) 

The function (|12f> reaches maximum at E x = 0.3736818 = Eco- In l° w fields 
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(\E X \ -C Eco) we have 



i. = (i - x [±^tim\ Ex = (pa\ Exi (13) 

J \ V 7tT erf(l/\/2T) J \T/ *' V J 

where angle brackets mean averaging over the equilibrium distribution. Note 
that the mobility temperature dependence in low fields (the expression within 
round brackets in Eq. is close to the analogous dependence for the mini- 

band cosine model (Ii(1/2T)/Iq(1/2T), I n (z) being the modified Bessel func- 
tion). 

Note, that the current maximum Ec(T) shifts leftwards (at E x > 0) and 
tends to its limiting value E\ — 0.2930031 with rising temperature. Essentially, 
the Ec value does not depend on the temperature at all in the cosine model. 
Meanwhile, in experiment [B] a shift of Ec to lower fields has been found with 
raising temperature. (In [B] miniband transport in 1SL GaAs/AlAs was in- 
vestigated in temperature range 10-300 K under fields up to 1.5 kV/cm, and 
negative differential conductivity was found surely enough.) 

We investigated numerically the temperature dependence Ec (T) . That func- 
tion can be approximated (with precision in calculations of 10~ 4 ) with the fol- 
lowing formula: 

E C {T) = \j E\ + (E ca - E\) cxp(-/3T) (E x - const, E C o = E(0)), (14) 

where (3 = 4.2. In Fig. the numerical solution (Ec(T)) of the equation 
dj x /dE x — is shown, the dotted curve being the approximation 114(1. 




Figure 3: Temperature dependence of the current maximum (Eq. I(lljl ). The 
dotted line shows approximation 1(14(1 . 
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3 Nonequilibrium phase transitions in 2SL 



The dispersion law of the square 2SL mentioned in Section Stakes the form 

p 2 + p 2 

g(p) = go + v , -irh/d < p X:V < nh/d, (15) 

where £q = const is electron energy in the lowest size-quantized level, A is 
miniband width. 

If the coordinate axes are directed along the SL principal ones, then the 
current density along X axis is given by Eq. I|ll|) . while that along Y axis is 
given by the same formula with x — > y change. 

Now a situation is considered, where the coordinate axes are rotated by angle 
45° with respect the principal axes of the square 2SL. In this case we have 



j y = E y + 4erf ( 1/V2T) sinh {1/(E X + E y )) 



erf 



erf 



exp 

VT 



T 



2(E X + E y y 
1 N 



erf 



\V2(E X + E y ) V2T 
4crf (l/Vzf) sinh (1/(E X - E y )) 

VT l N 

V2(E X -E y ) 



V2{E X +E y ) VW J 



erf 



exp 

VT 



T 



2(E X - Eyf 

1 \ 



V2(E X - Ey 



2T 



(16) 



where the field is expressed in units of V^^fi/edr, while the current in units of 
VZnh/edT. The formula for j x is obtained from Eq. i|16|) by interchange y «-> x. 
Let the sample is open in Y direction, so that 



jy = 0. 



(17) 



Then, by substituting Eq. I|16|) into Eq. (|17|) . we obtain an equation to find the 
spontaneous transverse field E y . At first, consider the case of extremely low 
temperatures (T — * 0) . In that case we get from Eq. I|16fl 



jy — E y ^ 



_sinh(l/(£; a; + Ey)) smh{l/{E x - E y )) 
The solution of Eq. (|17f) in reference to the current (|18fl has the form 



E u 



0, 
± 



< \EJ < Eco, 



coshq/^) 
sinh a (l/B a: ) 



E 2 C01 \E X \>E C0 



(18) 



(19) 



We write down, for comparison, the analogous formula for the 2SL cosine 
model: 

"0, < \E X \ < 1/tt, 



En, 



±VE 2 X - \E x \>l/n 



(20) 



6 



Figure 4: The transverse field E ys . The solid line shows numerical solution of 
Eqs. ftTjl. fT%| : the dotted one shows function fT^I . 

The numerical solution of Eq. (|17|l and function ljl9|l (the dotted curve) are 
shown in Fig. 0] 

The solution stability against small enough fluctuations of E y field is deter- 
mined with the following inequality [Sj: 



dj y _ <9 2 $ 



dE„ 



8E* 



>0, 



(21) 



where $ = J j y dE y + const is a synergetic potential (the entropy produc- 
tion PHI)- Conditions (fTTjl and (|2*T|l are satisfied at E y — E ys . Therefore, 
the solutions (|19|l and (|20|l correspond to minima of the potential $, which 
becomes double-well one (Fig. |SJ). Thus, there is a NPT2 here. 




: 5 



J i ---i i j_ 



4 p o 



Figure 5: Synergetic potential (up to a constant): a, - E x — Eco', b - E x = 1; c 
- = 1 for the cosine dispersion law. 



At T 0, the transverse field E y behaves the same way as at T — > with the 
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only difference that the bifurcation point Ec shifts leftwards with rising T (for 
E x > 0), in complete correspondence with behavior of the current j x maximum 
(see Fig. ©). 

It follows from the temperature dependence of Ec(T), that a novel type 
NPT2 is possible in the model considered. Indeed, in the field interval 0.317 < 
l-^zl < Eco the transverse field E ys (T) obeys the following relationship: 



E V8 (T) = ±J(E 2 CQ - El) [exp (-(3T C ) ~ exp (-0T)] (E x = const), (22) 



where the critical temperature is 



Tn = 4 In 



El 



P \cosh(l/E x )/ sixth 2 {1/E x )-Ef 



(23) 



Near the transition point Tc, the transverse field behaves as E ys 
The NPT2 is illustrated in Fig. El 
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Figure 6: The transverse field -B„ s at various values of the driving field: a 
E x = 0.318; b - E x = 0.330. 



4 Conclusion 

In present paper, an exact distribution function has been found of the carriers in 
the lowest parabolic miniband of 1SL. The novel formula for the current density 
in 1SL contains temperature dependence, which leads to the current maximum 
shift to the low field side with rising temperature, that agrees with available 
experimental data. 

The considered SL model includes NPT's studied earlier. Besides, a novel 
type NPT has been found, in which the sample temperature plays a role of a 
control parameter. 

The estimates of the effects predicted are reduced, in general, to estimate of 
electric field required. At d = 10~ 7 cm, r = 10~ 12 c, we have E a i=s 2 x 10 3 V/cm. 
At E x w 600 V/cm, the value T c « 45 K is obtained. 
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